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I. ABSTRACT 

We compute the probability of photon-axion oscilla- 
tions in the presence of both intergalactic magnetic fields 
and an electron plasma and investigate the effect on Type 
la supernovae observations. The conversion probability is 
calculated using a density matrix formalism by following 
light-paths through simulated universes in a Monte-Carlo 
fashion. We find, that even though the effect is highly 
frequency dependent, one needs to analyze relatively nar- 
row spectral features of high redshift objects in order to 
discern between the dimming effect from oscillations and 
a cosmological constant, in contrast to earlier claims that 
broad-band photometry is sufficient. 

PACS numbers: 13.40.Hq, 98.80.Es, 97.60.Bw 



II. INTRODUCTION 

Recently, Csaki, Kaloper and Terning Q| (CKT) pro- 
posed that the observed faintness of high redshift super- 
novae (SNe) could be attributed to the mixing of photons 
with a light axion in an intergalactic magnetic field. The 
nature of the oscillations is governed by the strength of 
the coupling which in turn depends on the axion coupling 
constant mass scale, the photon energy and the strength 
of the magnetic field. Assuming magnetic domains with 
uncorrelated field direction of size ~ Mpc, CKT found 
that for optical photons the oscillation is maximal and 
independent of energy, i.e., in the limit of infinite travel 
distance one approaches an equilibrium between the two 
photon polarization states and the axion. For optical 
photons, the probability to detect a photon as a function 
of the traveled distance, I, was approximated as 



■exp(-l/l ), 



(1) 



where Iq is the exponential decay-length. It was claimed 
that since this effect will cause additional dimming of 
high redshift SNe, constraints on the equation of state 
parameter of the dark energy component could be signif- 
icantly relaxed. (The presence of a non-clustering compo- 
nent of "dark energy" has recently been independently in- 
ferred from a combination of measurements of the cosmic 
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microwave background and the distribution of galaxies on 
large scales. The proposed photon-axion mixing does not 
remove the need for such a component, but its equation 
of state need not be as close to that for a cosmological 
constant as is the case without such a mixing.) For low 
photon energies, the oscillation was found to be energy 
dependent and the mixing very small. Therefore it should 
not severely affect the cosmic microwave background ra- 
diation which is redshifted to low energies at low z where 
the magnetic fields appear. For a magnetic field strength 
of \B\ ~ 1CP 9 G, CKT found that the current data can 
be accommodated by Hm = 0.3, fix = 0.7, wx = —1/3 if 
the axion mass is m ~ 10 -16 eV and the coupling scale 
is M ~ 4 ■ 10 11 GeV. 

It was then pointed out by Deffayet, Harari, Uzan and 
Zaldarriaga || that one has to take the effect of the in- 
tergalactic plasma into account, i.e., the free electrons, in 
which the photons are propagating. For a mean electron 
density of n c ~ 10 -7 cm -3 , they found that this will al- 
ter (lower) the oscillation probability and may also cause 
the oscillations for optical photons to be frequency de- 
pendent to a degree where the effect can be ruled out by 
observational constraints, namely by studying the color 
excess between the B and V wavelength bands. Only for 
very definitive properties of the intergalactic magnetic 
field could one get a large enough mixing angle and weak 
enough frequency dependence, namely with \B\ ~ 10~ 8 
G over domains of size ~ 10 kpc and weaker magnetic 
field strength over larger domains. (The strength and 
spatial properties of intergalactic magnetic are not well 
constrained by present measurements, for a recent review 
see §.) 

In a reply, Csaki et al. ||] pointed out that the mean 
electron density in most of space realistically is lower than 
the estimate used by Deffayet et al. by a factor of at least 
15. Since the energy dependence of the oscillations is 
very sensitive to this value, they find that this factor is 
enough to bring the energy dependence within current 
experimental bounds. 

In this note, we perform a full density matrix calcula- 
tion of the photon-axion oscillations and find the conver- 
sion probability to be highly frequency dependent and 
also not necessarily monotonically increasing with in- 
creasing redshift. We also calculate the color excess be- 
tween different wavelength bands for Type la SNe by inte- 
grating over spectrum templates modified by performing 
the density matrix calculation over the appropriate fre- 
quency range (taking the redshifting of the spectra into 
account). 
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III. DENSITY MATRIX FORMALISM 

We compute the mixing probability using the formal- 
ism of density matrices (see, e.g., 0). Following the no- 
tation of ||, we define the mixing matrix as 



M 




m sin a 



Am cos a 
Am sin a 
A m 



(2) 



The different quantities appearing in this matrix are 
given by 



A ± 
Ay =A ± , 

A M = 2 x 1CT 26 
A m = -2.5 x 10 



cm -1 , (3) 
(4) 



Bo,± \ ( Mg \ 1 _x (K , 



-28 



(™) 2 (^)"'™-'. < 6 » 



where 5o,j_ is the strength of the magnetic field perpen- 
dicular to the direction of the photon, M a is the inverse 
coupling between the photon and the axion, n c is the 
electron density, m a is the axion mass and u> is the en- 
ergy of the photon. The angle a is the angle between the 
(projected) magnetic field and the (arbitrary, but fixed) 
perpendicular polarization vector. Our standard set of 
parameter-values is given by 



Bo 



l(T 9 (l + z) 2 G, 
M a = 10 11 GeV, 
m a = l(T 16 eV, 
n c (z) = 1(T 8 (1 + z) 3 cm" 



(7) 



with a 20 % dispersion in B and n c and the redshift de- 
pendence of the magnetic field strength and the electron 
density comes from flux conservation and cosmological 
expansion, respectively. The equation to solve for the 
evolution of the density matrix p is given by 



with initial conditions 



Po 




(8) 



(9) 



Here the three diagonal elements refer to two different po- 
larization intensities and the axion intensity, respectively. 
We solve the system of 9 coupled (complex) differential 
equations numerically by following individual light 
paths through a large number of cells where the strength 
of the magnetic field and the electron density is deter- 
mined from predefined distributions and the direction of 



the magnetic field is random. Through each cell the back- 
ground cosmology and the wavelength of the photon are 
updated, as well as the matrices p and M. 

In order to study the qualitative behavior of the solu- 
tions, we rewrite M as a 2 x 2 matrix, 



M 2D = 



A A M 
Am A m 



(10) 



where A = Aj_ = A|| and Am is the component of the 
magnetic field parallel to the some average polarization 
vector of the photon beam. We solve Eq. || for the density 
matrix p 2D with initial conditions 



Po 



1 




(11) 



where the diagonal elements refer to the photon and the 
axion intensity respectively. Assuming a homogeneous 
magnetic field and electron density, we can solve the two- 
dimensional system analytically. For the p 2 ^ component, 
referring to the photon intensity, we get 



) 2 + 2A 2 



M 



(12) 



We see that we get maximal mixing if A = A, 



.... i-e., 

if m a = m max « 38 v /n o /(10- 8 cm- 3 ) 10- 16 eV. For 
m a 3> m max , the oscillations are suppressed as m a 4 . For 
m a <S TO max , the effect is insensitive to the values of the 
axion mass. Our numerical simulations show that even 



for m„ 



results are insensitive to the exact value 



of the axion mass. For values close to the typical set of 
parameter-values, we can set 



A 

27J' 



to get 



o 2D 
Pn 



A M 
A 



(1 



At 
COS 2^ ) - 



(13) 



(14) 



For small mixing angles, the effect should be roughly pro- 
portional to Bq and M~ 2 , whereas the effect should be 
rather insensitive to the exact values of the input parame- 
ters in cases of close to maximal mixing. We also expect 
the effect to be stronger for low values of the electron 
density approaching maximal mixing for ri = 0. These 
predictions are confirmed by numerical simulations. The 
oscillation length is of the order ~ Mpc. 



IV. RESULTS 

In Fig. [l], we show the attenuation due to photon-axion 
oscillations as a function of wavelength for one specific 
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line-of-sight for our standard set of input parameter- 
values (see Eq. 0) at z = 0.1 (upper panel), z = 0.2 
(middle panel) and z — 0.5 (lower panel) in a Om = 
0.28, Ha — 0.72 universe, which is what we will use sub- 
sequently. We have performed a number of simulations 
using a wide range of cosmological parameter-values and 
found the oscillation effect to depend only weakly on cos- 
mology. The most dramatic effect is the strong variation 
of attenuation with photon energy. 
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FIG. 1. The attenuation due to photon-axion oscillations 
for the standard set of parameter- values (see Eq. as a func- 
tion of wavelength at redshift z = 0.1 (upper panel), z = 0.2 
(middle panel) and z = 0.5 (lower panel). 

Since the attenuation varies very rapidly with photon 
energy in a similar manner over a broad energy range, 
we expect the frequency dependence to wash out to large 
extent when doing broad-band photometry. 

In Fig. ||, we show the rest-frame 5-band magnitude 
attenuation for Type la SNe due to photon-axion oscil- 
lations for three different values of the electron density, 
in the redshift interval < z < 2, using values for the 
other input parameters from Eq. |7|. Each point repre- 
sents the average value and the error bars the disper- 
sion for ten different lines-of-sight. In the upper panel, 



we have used n c — 
panel n c — 5 x 10~ 8 < 



10~ 7 cm- 
m~ 3 (l + 



! (1 + z) 3 , in the middle 
and in the lower panel 
n e = 10 _8 cm _a (l + z) d . Note that the effect is not nec- 
essarily increasing with increasing redshift. This is due 
to the fact that we are studying the rest-frame B-band 
magnitude. Since the amplitude of the oscillations scales 
roughly as (oj/n c ) 2 (see Eq. |lj), we need this combina- 
tion to be large at some point in order to get close to 
maximal mixing. If the plasma density is high, the pho- 



ton energy will be redshifted to too low energies before 
the plasma density is diluted due to the expansion. 
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FIG. 2. The attenuation integrated over the rest-frame 
B-band of Type la SNe due to photon-axion oscilla- 
tions with n — 10~ 7 cm~ 3 (l + z) 3 (upper panel), 
n c = 5 x 10" 8 cm" 3 (l + z) 3 (middle panel) and 
n c = 10 -8 cm~ 3 (l + z) 3 (lower panel). All other parame- 
ter-values are given by Eq. (?]. 

In Fig. ||, the rest-frame -B-band magnitude atten- 
uation for Type la SNe for three different values of 
the intergalactic field strength is shown. Again, each 
point represents the average value and the error bars 
the dispersion for ten different lines-of-sight. In the up- 
per panel, Bq = 10~ 10 (1 + z) 2 G, in the middle panel 
B = 5 x 10~ 10 (1 + z) 2 G and in the lower panel B = 



io- 9 

Eq 



i 1 



' G. All other parameter- values are given by 
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FIG. 3. The attenuation integrated over the rest-frame 
B-band of Type la SNe due to photon-axion oscillations with 
Bo = lCT 10 (l + ,z) 2 G (upper panel), B = 5 x l(T 10 (l + ,z) 2 G 
(middle panel) and B = 1CT 9 (1 + z) 2 G (lower panel). All 
other parameter-values are given by Eq. [?]. 

Our results indicate that we need low electron densi- 
ties in order to get an attenuation that increases with 
redshift. One should keep in mind that the overall nor- 
malization can be set by varying the strength of the mag- 
netic fields and/or the photon-axion coupling strength. 

Based on a sample of 36 z ~ 0.5 SNe Perlmutter et 
al. H measured the the average rest-frame B — V color 
excess to be 0.035 ± 0.022 mag. We find the expected 
color excess with our standard set of parameters to be, 
E(B — V) = 0.006, with a scatter around this value of 
0.004 mag, as shown in Fig. || where we have simulated 
200 Type la SNe at z=0.5. 
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FIG. 4. Color extinctions for rest-frame E(B - V) for 200 
Type la SNe at z = 0.5. 

We can investigate the frequency dependence when do- 
ing broad-band photometry by studying the extinction 
(at maximum intensity) in V — J, R — J and I — J for 
Type la SNe as a function of redshift, as we show in 
Fig. |[ All the broad-band filters and spectroscopy wave- 
length scales are in the observer's frame. We can see that 
in all three cases, the color excess is very small and thus 
difficult to measure. 
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FIG. 5. Color extinctions E(V — J), E(R-J) and E(I-J) 
for Type la SNe as a function of redshift. Note that the dim- 
ming associated with photon-axion mixing, unlike extinction 
by dust, also can generate blueing. 

In general, the nature of the frequency dependence will 
depend on the exact values of the input parameters with 
the possibility of generating both reddening and blueing, 
where in some cases the nature of the effect will depend 
on redshift (see, e.g., lower panel of Fig. ||). We thus 
conclude that it will be very difficult to use the color 
excess between different broad-bands to put severe limits 
on the photon-axion mixing parameters. 



ACKNOWLEDGEMENTS 

We are grateful to S. Hansen for bringing the CKT 
paper to our attention, and to G. Raffelt for insightful 
comments improving the quality of the manuscript. We 
also wish to thank J. Edsjo, H. Rubinstein, C. Fransson 
and J. Silk for helpful discussions. The research of L.B. is 
sponsored by the Swedish Research Council (VR). A.G. is 
a Royal Swedish Academy Research Fellow supported by 
a grant from the Knut and Alice Wallenberg Foundation. 



V. DISCUSSION 



Our numerical simulations indicate that in order to 
get a dimming effect from photon-axion oscillations sim- 
ilar to the one from a cosmological constant (increas- 
ing at lower redshifts, saturating at higher), one would 
need to have an average intergalactic electron density of 
n e ^5 10~ 8 cm~ 3 (l + z) 3 . Assuming this, it should be 
possible to vary the average magnetic field strength or 
the photon-axion coupling strength in order to fit the 
current broad-band photometry data. Note that in the 
case of close to maximal mixing, results are generally not 
very sensitive to the exact values of the input parameters, 
yielding results similar to the upper panel in Fig. |^. 

Depending on the choice of mixing parameters the 
color excess terms of sources at cosmological distances 
could be either positive or negative, the latter case being 
particularly interesting as it could not be confused with 
regular extinction by dust. Since photon-axion oscilla- 
tions can cause either reddening or blueing (or no color 
excess at all) for close to maximal mixing and the inte- 
grated broad-band magnitudes wash out the dispersion 
in attenuation, we expect spectroscopic studies of high- 
z objects to be a more powerful discriminator between 
different oscillation models and the case of a cosmolog- 
ical constant or any other dark energy component (for 
which we suppose the dimming to be entirely frequency 
independent). Systematic analysis of quasar, gamma-ray 
burst and galaxy spectra as a function of redshift by, 
e.g., the Sloan Digital Sky Survey and 2dF groups are 
probably the best probes for the photon-axion mixing 
parameter space. Note that the source size sets a lower 
limit to the size of the fluctuations that can be probed 
since the fluctuations will average out if photons from 
different parts of the source will travel through different 
magnetic field strengths and electron densities. 

SN spectra, available with good signal-to-noise ratios 
up to at least z ~ 0.5, are useful probes for the large 
mixing parameters that would yield a sufficient dimming 
of Type la SNe as to explain the Hubble diagrams in [|]|]] 
without invoking dark energy. 
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